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The charged pion form factor, FTy{Q ), is an important quantity which can be used to advance our 
knowledge of hadronic structure. However, the extraction of from data requires a model of the 
^E[(e, e'7r+)n reaction, and thus is inherently model dependent. Therefore, a detailed description 
of the extraction of the charged pion form factor from electroproduction data obtained recently 
at Jefferson Lab is presented, with particular focus given to the dominant uncertainties in this 
procedure. Resuhs for F.^ are presented for Q^=0. 60-2.45 GeV^. Above = 1.5 GeV^, the F^ 
values are systematically below the monopole parameterization that describes the low data used 
to determine the pion charge radius. The pion form factor can be calculated in a wide variety 
of theoretical approaches, and the experimental results are compared to a number of calculations. 
This comparison is helpful in understanding the role of soft versus hard contributions to hadronic 
structure in the intermediate regime. 

PACS numbers: 14.40.Aq,13.40.Gp,13.60.Le,25.30.Rw,11.55. Jy 



I. INTRODUCTION of their constituents, the quarks and gluons. However, 

this structure is too complicated to be calculated rig- 
There is much interest in trying to understand the 
structure of hadrons, both mesons and baryons, in terms 
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orously in Quantum Chromodynamics (QCD) because 
perturbative QCD (pQCD) methods are not applicable 
in the confinement regime. Chiral Perturbation Theory 
can give valuable insights, but it is limited to small val- 
ues of the photon virtuality . Hence, in the inter- 
mediate regime one has to resort to models like the 
constituent quark model or methods employing Light- 
Cone (LC) dynamics or the Bethe-Salpeter (plus Dyson- 
Schwinger) equation, or to other approaches such as the 
use of dispersion relations or (QCD or LC) sum rules. 

Transitions and (transition) form factors are crucial 
elements for gauging the ideas underlying these QCD- 
based models. For example, the constituent quark model 
gives a fairly good description of the meson and baryon 
spectrum and some transitions, but quark effective form 
factors are typically required when describing hadronic 



form factors in the experimentally accessible region. 
In this framework, the study of hadronic form factors 
can thus be viewed as a study of the transition from 
constituent to current quark degrees of freedom. As ex- 
emplified by the many calculations of it, the electric form 
factor of the pion, i^^r, is one of the best observables for 
the investigation of the transition of QCD effective de- 
grees of freedom in the soft regime, governed by all kinds 
of quark-gluon correlations at low Q^, to the perturba- 
tive (including next-to-leading order and transverse cor- 
rections) regime at higher . 

In contrast to the nucleon, the asymptotic normaliza- 
tion of the pion wave function is known from pion decay. 
The hard part of the 7r+ form factor can be calculated 
within the framework of pQCD as the sum of logarithms 
and powers of [1] 



i^.(Q') 



47rCFa,(Q2) 



r 



[l + 0(a,(Q2),^/Q2)]^ 



(1) 



which in the ^ oo limit becomes [l|, 



(2) 



where f-^ = 93 MeV is the pion decay constant Q . 

Because the pion's qq valence structure is relatively 
simple, the transition from "soft" (non-perturbative) to 
"hard" (perturbative) QCD is expected to occur at sig- 
nificantly lower values of for Fj, than for the nucleon 
form factors Q. Some estimates [d] suggest that pQCD 
contributions to the pion form factor are already signifi- 
cant at > 5 GeV^ . On the other hand, a recent analy- 
sis [6*1 indicates that non-perturbative contributions dom- 
inate the pion form factor up to relatively large values of 
Q^, giving more than half of the pion form factor up to 
Q^=20 GeV^. Thus, there is an ongoing theoretical de- 
bate on the interplay of these hard and soft components 
at intermediate Q^, and high quality experimental data 
are needed to help guide this discussion. 

In this work, we concentrate exclusively on the space- 
like region of the pion form factor. For recent measure- 
ments in the timelike region see Ref. Q • At low values of 
Q^, where it is governed by the charge radius of the pion, 
has been determined up to (3^=0.253 GeV^ %^ from 
the scattering of high-energy pions by atomic electrons. 
For the determination of the pion form factor at higher 
values of Q^, one has to use high-energy electroproduc- 
tion of pions on a nucleon, i.e., employ the ^II(e, e'7r+)?i 
reaction. For selected kinematic conditions, the longi- 
tudinal cross section is very sensitive to the pion form 
factor. In this way, data for F^^ have been obtained for 
values of up to 10 GeV^ at Cornell How- 
ever, those data suffer from relatively large statistical 



and systematic uncertainties. More precise data were ob- 
tained at the Deutsches Elektronen-Synchrotron (DESY) 
jl3, 14]. With the availability of high-intensity electron 
beams, combined with accurate magnetic spectrometers 
at the Thomas Jefferson National Accelerator Facility 
(JLab), it has been possible to determine L/T separated 
cross sections with high precision. The measurement of 
these cross sections in the regime of (5^=0.60-1.60 GeV^ 
[Experiment Fpi-l_[ll,[ll]] and Q^=1.60-2A5 GeV^ [Ex- 
periment Fpi- 2 1 1 Til l are described in detail in the pre- 
ceding paper [l8|. In this paper, it is discussed how to 
determine from measured longitudinal cross sections, 
the values determined from the JLab and DESY data are 
presented, and the results of various theoretical calcula- 
tions are compared with the experimental data. 



Since the pion in the proton is virtual (off its mass- 
shell), the extraction of F-^ from the measured electropro- 
duction cross sections requires some model or procedure. 
In the next section, the methods that have been used to 
determine F^^ from the data are discussed. Section IIIII 
presents the adopted extraction method and the values 
of F^ thus determined, including a full discussion of the 
uncertainties resulting from the experimental data and 
those from the adopted extraction procedure. Various 
model calculations of are discussed and compared to 
the data in section HVl In the final section, some conclu- 
sions are drawn and an outlook for the future is given. 
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II. METHODS OF DETERMINING THE PION 
CHARGE FORM FACTOR FROM DATA 

The measurement of the pion form factor is chaUeng- 
ing. As stated in the introduction, at low can be 

measured in a model-independent manner via the elastic 
scattering of tt"*" from atomic electrons, such as has been 
done up to Q^=0.253 GeV^ at Fermilab [8] and at the 
CERN SPS It is not possible to access significantly 
higher values of with this technique because of limi- 
tations in the energy of the pion beam together with the 
unfavorable momentum transfer. Therefore, at higher 
values of F-^ must be determined from pion electro- 
production on the proton. The dependence on F^. enters 
the cross section via the ^-channel process, in which the 
incident electron scatters from a virtual pion, bringing it 
on-shcll. This process dominates near the pion-pole at 
t — m'^, where t is the Mandelstam variable. The phys- 
ical region for t in pion electroproduction is negative, so 
measurements should be performed at the smallest at- 
tainable values of —t. To minimize background contribu- 
tions, it is also necessary to separate out the longitudinal 
cross section ctl, via a Rosenbluth L/T(/LT/TT) sepa- 
ration [T9| . 

The minimum physical value of —t, —tmim is non-zero 
and increases with increasing and decreasing value of 
the invariant mass, W, of the produced pion-nucleon sys- 
tem. Carlson and Milana [20| have estimated an approxi- 
mate upper limit for the value of —tmin of the data appro- 
priate for the extraction of the pion form factor by study- 
ing the competing non-pole QCD processes, which may 
complicate the extraction of at higher . They found 
that the background ratio MpQco/Mpoie rises dramati- 
cally once —t„iin > 0.20 GeV^. Their concern stemmed 
from the large value of —t in some of the Cornell results, 
which have —tmin > 0.4 GeV^ Therefore, reliable 

Ftt measurements should be performed at smaller —t and 
thus higher W (for a fixed Q^). The results presented in 
this paper respect this —tmin < 0.20 GeV^ upper limit. 
It is yet to be determined if reliable F-^ measurements 
can be made in the future at larger —t. 

The value of F^(Q^) can then be determined from the 
data by trying to extrapolate the measured longitudi- 
nal cross sections at small values of —t to the pole at 
t = = 0.02 GeV^, or by comparing the measured lon- 
gitudinal cross section at small values of —t to the best 
available model for the ^H(e, e'7r"'")n reaction, adjusting 
the value of in the latter. The presence of the nucleon 
and its structure complicates the theoretical model used, 
and so an unavoidable implication of this method is that 
the extracted pion form factor values are model depen- 
dent. The differential cross sections ctl versus t over some 
range of and W are the actual observables measured 
by the experiment. It is important to note that in all 
cases the use of a model to extract is justified only if 
the model correctly predicts the t-dependence and mag- 
nitude of the (Tl data as well as the dependence on the 
invariant mass W of the photon-nucleon system. 



A. Chew-Low Extrapolation Method 

Frazer [2l[ originally proposed that F^. be extracted 
from (Tl via a kinematic extrapolation to the pion-pole, 
and that this be done in an analytical manner using the so 
called Chew-Low extrapolation [S^l ■ The used Born for- 
mula is not gauge invariant p3| . but in principle should 
give Fn-, nonetheless, when extrapolating to the pole. 

The last serious attempt to extract the space-like pion 
form factor from electroproduction data via the Chew- 
Low method was by Devenish and Lyth in 1972. 
Most of the data used were unseparated cross sections. 
The extrapolation failed to produce a reliable result, be- 
cause different polynomial fits that were equally likely 
in the physical region gave divergent values of the form 
factor when extrapolated to the pion-pole at i = m^. 
Since then, the quality of the tt"*" electroproduction data- 
set has improved immensely, and separated longitudinal 
cross sections can now be used, avoiding the complica- 
tions stemming from the other parts of the cross section. 
Therefore, it has been suggested to us that it may be ap- 
propriate to revisit the Chew-Low extrapolation method. 

However, before trying this method on the new data, it 
should be tested to see how reliably one can extrapolate 
to the pole. We start with high precision ctl 'pseudo- 
data' generated as a function of —t with the VGL Regge 
model. This model gives a fair to good description of 
a wide body of pion photo- and electroproduction data 
(see section Hi CI) . The kinematic conditions for the test 
are ^ 1,594 Q^y^ and = 2.213 GeV, similar to 
our Fpi-2 data. The input value of the pion form factor 
in the model was F^r = 0.244. The model ctl cross sec- 
tions were then used in a Chew-Low type extrapolation, 
with the challenge being to see if the Chew-Low extrap- 
olation is able to reproduce (within fitting uncertainties) 
the input F^- -value. 

The basis of the Chew-Low method is the Born-term 
model (BTM) formula for the pion-pole contribution to 
(Tl. We use the BTM of Actor, Korner and Bender [2^, 
where pion-pole contribution to ctl is given by 

= 4hc {eg^NN? J^^^ Q'F^{Q'), (3) 

where e^/(47r?ic) = 1/137 and N, which depends on the 
flux factor used in the definition of is given in our 
case by 

N = 32n(W^-ml)^(W^ - + + 2Q'^{W^ + m^) 

. . 

|26l |27[ . A monopole parameterization of the g-^NN form 
factor is typically used to determine its value at i-values 
away from the pion-pole 

where g-wNNirn^) is the experimental value of 13.4 [2^ . 
This is also the value used in the VGL calculations. We 
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FIG. 1: (Color online) Linear (dotted), quadratic (dashed) 
and cubic (solid line) extrapolations of to the pole as com- 
puted from Eqn. [6l The boxes are a VGL Regge model calcu- 
lation for (JL at fixed W = 2.213 GeV and = 1.594 GeV^ 
calculated with ,^^=0.244. The lower limit of the box range 
is the kinematic endpoint of these Q^, W values, while the 
upper limit is given by the f-range of our experiment. The 
input Ftt value in the model is indicated by the bullet placed 
at the pion-pole. 



use the = 0.80 GeV result from the QCD Sum Rules 
calculation by T. Meissner [29|], but because of the ex- 
trapolation to the pole the final result does not depend 
significantly upon the value chosen. 

For the Chew-Low extrapolation, one plots the value 

of 

j^2_ N (t - mlf daL 

versus ~t, which for a pure pole cross section gives 
a straight line passing through the origin, with value 
Ft^IQ"^) at the pole it = mj). Other contributions to the 
cross section, which have to be present, because the pole 
contribution alone is not gauge invariant, will change this 
behavior, but since they do not contain the fac- 
tor, they will not influence the value of at the pole. 
However, it is not a priori given that the behavior as 
function of —t is linear, quadratic, or of higher order, 
thus introducing a 'model' (extrapolation) uncertainty. 



Values of for the generated pseudodata, together 
with linear, quadratic and cubic extrapolations to the 
pole are shown in Fig. [TJ Also shown is the input form 
factor value in the VGL model, plotted at the pion-pole. 
Quadratic and higher-order extrapolations are almost in- 
distinguishable and give a very good description of the 
(pseudo)data, but miss the input value of F-^. This was 
true for all cases that were investigated, from = 0.60 
to 2.45 GeV^, the deviation from the input Fir-value be- 
ing 6-15%, depending on the case and the order of the 
extrapolation polynomial. Overall, there was no consis- 
tent trend for the order of polynomial which was best 
able to reproduce the input form factor value. 

This study indicates that even if ctl is very well known 
over a range of physically-accessible t, the Chew Low ex- 
trapolation yields inconsistent results. The extrapolated 
result depends greatly upon the choice of quadratic cubic, 
or higher-order function, which all give a very good de- 
scription of the data in the physical region. This indicates 
that the t-dependence of data in the physical region is in- 
sufficient to uniquely constrain the extrapolation through 
the unphysical region to the pole, even if the data have 
small relative uncertainties. Furthermore, even though 
modern data such as the JLab (Tl data are much more 
precise than those previously available, they still com- 
prise 4-6 i-bins only, each with statistical and systematic 
uncertainties of 5-10%. Therefore, any polynomial ex- 
trapolation of such data to the pole will be more unreli- 
able than the pseudodata test case shown here. There- 
fore, the Chew-Low extrapolation technique cannot be 
used to reliably determine the pion form factor from a 
realistic ctl data set. 



B. Early Extractions of F,r 

Brown et al. 30] at CEA were the first to embrace the 
use of theoretical input to determine F^- from their data. 
They used the model of Berends |31| , which includes the 
dominant isovector Born term, with corrections for t val- 
ues away from the pole by means of fixed-^ dispersion 
relations. This model was also used by Bebek et al. for 
the analysis of the first two sets of Cornell data [13, [HI ■ 
The model gave a fair description of the data, but system- 
atically underpredicted the LT term of the cross section 
and the i-dependence of the data. 

Until then, data were obtained at one (larger) value 
of the photon polarization parameter e only. In the 
third Cornell experiment [12], data were taken at low 
values of e, so that in combination with the earlier data 
an L/T separation could be performed at Q^-values of 
1.19, 2.00 and 3.32 GeV^. The value of ctt was found 
to be substantially larger than predicted by Berends, es- 
pecially at larger Q"^ . The values obtained for ctl had 
such large error bars that they were not used to deter- 
mine Fn-. Instead, use was made of the observation that 
within the experimental error bars the Q^-dependence 
of the forward transverse cross section was satisfactorily 
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reproduced by the Q^-dependence of the total virtual- 
photoproduction cross section. Therefore, crT(Q^) was 
parameterized with the overall scale as a free parameter, 
and the parameterized values then used to subtract ctt 
from the measured unseparated cross sections to obtain 
ctl- These ctl data at the lowest value of —t were used 
to determine F-^ , assuming that ctl is given there by the 
i-channel one pion-exchange Born term. This was done 
for all data obtained at CEA and Cornell. The uncer- 
tainties in thus obtained and presented in Ref. [l^ 
are statistical ones only, and do not include the contribu- 
tion from the uncertainty in the value of ctt used in the 
subtraction. Especially at the larger values of Q^, these 
are considered to be substantial, as can be seen from Fig. 
4 of Ref. 0. 

The DESY experiments produced high-quality sepa- 
rated cross sections at = 0.35 GeV^, W = 2.10 GeV 
[i3 and g2 = 0.70 GeV^, W = 2.19 GeV [3]. Both of 
these experiments used the generalized Born Term Model 
of Gutbrod and Kramer 3^ to determine Ft^. This 
BTM incorporates t, s, and u-channel diagrams for the 
7^ -l- p — > TT"*" -|- n reaction, giving a fair description of the 
magnitude of the measured unseparated cross sections, 
but failing to describe cttt and ctlt- However, Gutbrod 
and Kramer found that when treating the magnitude of 
the nucleon form factor G^{Q^) as a free parameter, a 
much better description of the then available data was 
obtained. In addition, they included a factor e^^^'^ in 
order to improve the description of the t dependence of 
the data. The justification given is that the nucleon is 
far off its mass-shell, whereas the pion is near to its pole. 
This generalized BTM gave a good overall description 
of the DESY data. However, at = 0.70 GeV^, nu- 
cleon form factors about 50% above their on-mass-shell 
values were needed. The size of the modification needed 
at Q2 ^ 0.35 GeV^ is not given. 



C. Newer Models 

More recently, two new models for the ^H(e,e'7r+)n 
reaction have become available. 

In Refs. [H, 111], Vanderhaeghen, Guidal and Laget 
(VGL) have presented a Regge model for pion production 
in which the pole-like propagators of Born term models 
are replaced with Regge propagators, i.e., the interaction 
is effectively described by the exchange of a family of par- 
ticles with the same quantum numbers instead of a single 
particle. If the same vertices and coupling constants are 
used, the Regge model and the BTM calculations agree 
at the pole of the exchanged particle, but away from 
the pole the Regge model provides a superior description 
of the available data. For forward pion production, the 
dominant exchanges are the tt and p trajectories. These 
determine the t-dependence of the cross section without 
the use of a g-nNNit) factor. At low values of —t, as cov- 
ered by this work, ctl is completely determined by the 
TT trajectory, while ut is also sensitive to the p exchange 



contribution. Since the t-channel tt diagram is by itself 
not gauge invariant, the s-channel (for 7r+ production) or 
M-channel (for tt^ production) nucleon exchange diagram 
was also Reggeizcd, to ensure gauge invariance of their 
sum. 

The VGL model was first applied to pion photopro- 
duction 33] and later extended to electroproduction [23[ , 
with monopole forms for the 77777 and ^777 form factors: 



i^..p(Q') = [i + OVA^J- 



(7) 



Apart from the 77777 and ^777 form factors, the model is 
parameter free, as the coupling constants at the vertices 
(such as Qp-n-y) are well determined by precise studies and 
analyses in the resonance region. The model gives a good 
description of the W- and ^-dependences of then avail- 
able 77+ and 77~ photoproduction data, including the spin 
asymmetries, and of the earlier electroproduction data. 

The VGL predictions have been compared to our mea- 
sured cross sections and the ones taken at DESY [l^ [l3| 
in Ref. ^IS^. For the discussion in this paper, the data 
for (Jl and (Tt are reproduced in Fig. [21 together with 
the results of the model calculations. The VGL cross 
sections were evaluated at the same W and Q values 
as the data, resulting in the discontinuities shown. The 
values of Aj shown are determined by the fitting of the 
VGL model to the measured ctl -values at the five values 
oit at each Q^, resulting in values between 0.37 and 0.51 
GeV^. The value of A^ is more poorly known. Galcu- 
lations with both A^=0.600 and 1.500 GeV^ are shown, 
where the upper value is taken from the application of 
the VGL model to kaon electroproduction [34| . 

The model gives an overall good description of our ctl 
data and those of [l^ but the description of the 
t-dependence of the data is worse at Q^=0.60 and 0.70 
GeV^. The poorer description of the ctl data by the VGL 
model at lower and W may be due to contributions 
from resonances, which are not included explicitly in the 
Regge model. This is supported by the fact that the dis- 
crepancy in the i-dependence of the ctl data is strongest 
at the lowest value, at higher the resonance form 
factor supposedly reducing such contributions. The val- 
ues of (Tt are severely underestimated, especially at larger 

, even when taking a hard /0777 form factor. Since 
the data at the real-photon point are well described, this 
suggests that another mechanism, whose contribution in- 
creases with Q\ is at play [35|] . Recently the VGL model 
was extended [3^ by including a hard scattering between 
the virtual photon and a quark, the latter hadronizing in 
combination with the spectator diquark into a pion plus 
residual nucleon. With plausible assumptions, a good 
description of ctt was obtained, with no influence on (Tl. 
Those results support the idea that the discrepancy in 
the magnitude of ctt, which increases with Q^, and the 
discrepancy in the slope of ctl with —t, which decreases 
with Q^, are not directly related. Strategies for deal- 
ing with the latter discrepancy when extracting the pion 
form factor are discussed in Sec. IIIII 
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Q'=0.60GeV2 W=1.95GeV Q^=0.75GeV2 W=1.95GeV Q^=1.00GeV^ W=1.95GeV q^=1,60GeV^ W = 1.95GeV 




-t (GeV^) 

FIG. 2: (Color online) Separated tt^ electroproduction cross sections ctl [solid] and ctt [open] from JLab and DESY in 
comparison to the predictions of the VGL Regge model [23l ]. The error bars of the JLab data represent the combination of 
statistical and t uncorrelated systematic uncertainties. In addition, there is an overall systematic uncertainty of about 6%, 
mainly from the t correlated, e uncorrelated systematic uncertainty. The VGL Regge model calculations for Q^=0.60-1.60 
GeV^ GeV use A^=0.394, 0.372, 0.411, 0.455, GeV^ and those for Q2=0.35-2.45 GeV^ W ~ 2.1 GeV use A^=0.601, 

0.519, 0.513, 0.491 GeV^ The solid(dashed) curves indicate the A^ = l. 500(0. 600) GeV^ value used. 
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We also considered a modification to the VGL Regge 
model published by J.M. Laget in 2004 |37|]. Laget in- 
troduces a t-dependent factor into the pion form factor 
which is related to the pion saturating Regge trajectory, 
approaching -1 as t — > — oo. The effect of this modifi- 
cation is to boost (Tt by 40% for the largest —t spanned 
by our data (Q2^2.703 GcV^, -i=0.365 GeV^), and con- 
verging with the unmodified calculation at small —t. The 
effect on (Tl is under 1% for the largest — t covered by our 
data, and is negligible at —tmin- 

Another recent development is the effective La- 
grangian model of Faessler Gutsche, Lyubovitskij and 
Obukhovsky (FGLO, Ref. |i, H^). This is a modified 
Born Term Model, in which an effective Lagrangian is 
used to describe nucleon, pion, p and photon degrees of 
freedom. The (combined) effect of s- and w-channel con- 
tributions, which interferes with the pion t-pole, is mod- 
eled using a constituent quark model. The authors show 
that the p i-pole contribution is very important in the 
description of the magnitude of ut. When comparing 
vector and tensor representations of the p contribution, 
the latter was found to give better results. Unlike the 
VGL model, the ctl cross section depends here also on 
the p exchange, because of the interference of the tt and 
tensor p exchange contributions. The model contains a 
few free parameters, such as the renormalization con- 
stant of the KroU-Ruderman contact term used to model 
the s(M)-channel, and t-dependent strong meson- nucleon 
vertices, which are parameterized in monopole form, as 
are the electromagnetic form factors. The correspond- 
ing parameters were adjusted so as to give overall good 
agreement with our ctl and ctt data. 

As in case of the VGL model, a detailed comparison of 
the FGLO model results to the measured data is given in 
Ref. ^IS*!, while the results for ctl and ctt are also shown 
in Fig. [3l The values of A^ used were determined by the 
fitting of the model to the ctl t-bins at each Q^, while 
keeping the other parameters fixed at the values assigned 
by the authors. In some cases^this results in different A^ 
values than shown in Ref. [33|. However, it should be 
kept in mind that the FGLO model ul cross sections also 
depend on other parameters, which have been adjusted 
by the authors of the model to give good agreement to 
our ctl and ctt data. To the best of our knowledge, the 
Q^=0.7 GeV^ data of Ref. [3| were not taken into ac- 
count when these parameters were determined. 

Generally, the agreement of the FGLO model with the 
(Tl data is rather good except for the (5^=0.60 [Fpi-1] and 
0.70 GeV^ measurements. There is a serious discrep- 
ancy in the Q^- and VF-dependence of the ctt data. For 
around 0.7 GeV^, the model agrees fairly well with the 
data atW = 1.95 GeV, but it over-predicts the (5^=0.70 
GeV^, W = 2.19 GeV data by a large factor. On the 
other hand, for Q^^LGO GeV^, the W = 1.95 GeV data 
are under-predicted by about a factor of two, while those 
at Vl^ = 2.22 GeV are reproduced, and the W = 2.22 
GeV data for Q^—2A5 GeV^ are under-predicted again 
by 20-60%. This indicates some problem in the descrip- 



0' w 

(GeV^) (GeV) 


A" 
(GeV^) 




0.60 1.95 
0.75 1.95 
1.00 1.95 
1.60 1.95 


0.458 ± Omito lll 
0.388 ± 0.038lJ5:o53 
0.454 ± 0.0341°:°;^^ 
0.485 ± 0.038io:o27 


0.433 ± O.Ontolll 
0.341 ± 0.022j:o.!J3? 

0.312 ± o.oielE! j;?^ 

0.233 ± 0.014t;j;^ig 


0.35 2.10 
0.70 2.19 


0.601 ± 0.060 
0.627 ± 0.058io:o85 


0.632 ± 0.023 
0.473 ± 0.023to;o3i 


1.60 2.22 
2.45 2.22 


0.513 ± 0.033lJ5:^g 
0.491 ± 0.035lJ5:o24 


0.243 ± 0.012tE;;°J^ 
0.167 ±0.010lE!;JjJ^ 



TABLE I: and values from this work, and the reana- 
lyzed data from Refs. [l3 . [Til ] using the same method. The 
first error includes all experimental and analysis uncertainties, 
and the second error is the 'model uncertainty' as described 
in the text. The total uncertainty is found by taking their 
sum, in quadrature. Please note that in some cases the A^ 
value listed is different than the value used in Fig. [21 

tion of the Q^, 14^-dependences of the p exchange used 
to describe ctt. Because of the p — tt interference, the 
problems with the description of ax also affect the ctl 
calculation. This makes it hard to estimate how reliable 
the values of would be if extracted from the data using 
this model. 



III. RESULTS 

As already discussed, the separated cross sections ver- 
sus t over some range of and W are the actual observ- 
ables measured by the experiment, and the extraction of 
the pion form factor from these data is inherently model 
dependent. Ideally, one would like to have a variety of 
reliable electroproduction models to choose from, so that 
the model dependence of the extracted F.^ values can be 
better understood. Since the VGL Regge model is able, 
without fitted parameters, to provide a good description 
of both TT"*" and tt~ photoproduction data, and of ctl elec- 
troproduction data over a range in W, t, and it IS our 
opinion that at the moment only this model has shown 
itself to be sufficiently reliable to enable its use to extract 
pion form factor values from the ctl data. Therefore, we 
will use this model to determine our F^^ values. Clearly, 
the i^TT values determined are strictly within the context 
of the VGL Regge model, and other values may result if 
other, better models become available in the future. 



A. W ^ 2.2 GeV Data 

As shown in Fig. [21 the VGL model does a good job 
of describing the i-dependence of the ctl cross sections 
at W « 2.2 GeV, =0.35, 1.60 and 2.45 GeV^. In 
these cases, the extraction of the pion form factor from 
the data is straightforward: the value of A^ in the model 
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FIG. 3: (Color online) Separated tt"*" electroproduction cross sections ctl [solid], and ctt [open] from this work and DESY [ij] 
in comparison to the FGLO effective Lagrangian model 39]. The data error bars and systematic uncertainties are as in Fig. 
El The solid (dashed) curves denote model calculations for (jl ((Jt) with A^=0.405, 0.414, 0.503, 0.654, 0.386, 0.608 and 0.636 
GeV^ (from upper left to lower right). The calculations were performed at the same W and Q as the data, with straight lines 
connecting the calculated values. 
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is varied until the agreement of the model with the data 
is optimized. The mean Q and W values of the data 
for each i-bin are used when evaluating the model. 
is then calculated from Eqn. [3 using the best-fit and 
the nominal values. These arc listed in the last two 
lines of Table m 

The experimental statistical and systematic uncertain- 
ties were propagated to the Ft, uncertainties as follows. 
The statistical and t, e-uncorrelated systematic uncer- 
tainties^ were applied to the ctl data prior to the fitting of 
the VGL model to the ctl data. This yields the best-fit 
value and its associated fitting uncertainty. The effects 
of the i-correlated, e-uncorrelated, and the t, e-correlated 
systematic uncertainties on the fit were determined by in- 
vestigating the variation in A^ values allowed by fitting 
to the lowest —t bin only. Of these, the e-uncorrelated, t- 
correlated systematic uncertainty is amplified by 1/ Ae in 
the L-T separation, while the t, e-correlated uncertainty 
is not. The resulting uncertainties are added in quadra- 
ture to the fitting error, yielding the first Aj uncertainty 
listed in Table HI This value is also propagated to Ft, ac- 
cording to the monopole parameterization, yielding the 
first Ft; uncertainty listed. 

In order to check if the extracted value of Ft, depends 
on the t-range used, the VGL model (i.e., the value of 
A^) was fitted separately to each ctl point from Fpi-2 
and DESY [l^, , and the corresponding values of Ft, 
determined. In order to remove the natural variation of 
Ft, with the Q of each bin, the nominal values were 
used in the monopole equation. A plot of the obtained Ft, 
versus t is shown in Fig. IH Also indicated as the shaded 
band is the Ft, value with the uncertainty that is obtained 
if one fits to all of the t-bins simultaneously. Except 
perhaps at (5^=0.70 GeV^, the data show no residual 
i-dependence beyond the statistical fluctuation. 
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FIG. 4: (Color online) Ft, consistency check for the DESY and 
Fpi-2 data at W ~ 2.2 GeV. The solid squares indicate the 
Ft, values that would be obtained if the VGL model was fit 
to each ctl point separately. The shaded band is the Ft, value 
that is obtained if the model is fit to all of the t-bins. The 
error bars and band refiect the statistical and t-uncorrelated 
systematic uncertainties only. 



B. W = 1.95 GeV Data 

As already shown in Sec. Ill CI the VGL model does not 
fully describe the ^-dependence of our ctl data at W = 
1.95 GeV. The difficulty, as far as the Ft, extraction is 
concerned, is that there is no theoretical guidance for the 
assumed interfering background not included in the VGL 
model, even if one assumes that it is due to resonances. 
Virtually nothing is known about the L/T character of 
resonances a.t W = 1.95 GeV, let alone how they may 
influence ctl through their interference with the 7r-pole 
amplitude. Given this lack of theoretical guidance, we are 
forced to make some assumptions in extracting Ft, from 
these data. Our guiding principle is to minimize these 
assumptions to the greatest extent possible. The form 
factor extraction method that we have adopted for these 
data relies on the single assumption that the contribution 



^ These uncertainties are described in detail in Ref. Il8|l 



of the background is smallest at the kinematic endpoint 

Our best estimate of Ft, for the W — 1.95 GeV data is 
determined in the following manner. Using the value of 
Aj as a free parameter, the VGL model was fitted to each 
t-bin separately, yielding A'^{Q'^, W, t) values as shown in 
Fig. [5l The values of A^ tend to decrease as —t increases, 
presumably because of an interfering background not in- 
cluded in the VGL model. Since the pole cross section 
containing Ft, increases strongly with decreasing — we 
assume that the effect of this background will be small- 
est at the lowest value of |f| allowed by the experimental 
kinematics, |tmm|. Thus, an extrapolation of Aj to this 
physical limit is used to obtain our best estimate of Ft,. 
The value of A^ at tmin is obtained by a linear fit to the 
data in Fig. O The resulting Aj and Ft, values for the 
Fpi-I data are listed in Table HI The first uncertainty 
listed includes both the experimental and the linear fit 
extrapolation uncertainties. 

Since Fig. [5] suggests also a dependence (at larger —t) 
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FIG. 5: (Color online) Values of determined from the fit of 
the VGL model to each t-hin, and linear fit to same. The error 
bars refiect the statistical and t-uncorrelated systematic un- 
certainties. The additional overall systematic uncertainties, 
which were applied after the fit, are not shown. 



between the VGL calculation and the Q^=0.70 GeV^ 
data of Ref. [14], this extraction method was also 
applied to those data. The result obtained when extrap- 
olating to tmin is hstcd in Table U The value of F7r(Aj) 
is 11(20)% larger than if the VGL model was simply fit 
to all data points. Applying the same procedure to our 
W = 2.22 GeV data, it was found that the resulting 
values of F7r(A^) would be 1(2)% larger, which is statis- 
tically insignificant, confirming that the i-dependence of 
those data is well described by the VGL model. 



C. Model Uncertainty Estimate 

The fact that we used an additional assumption for 
the cases where the VGL model does not completely de- 



scribe the i-dependence of the (Tl data causes an addi- 
tional uncertainty in the extracted F^^ value, which we 
term 'model uncertainty'. This model uncertainty, which 
is within the context of the VGL model, should be dis- 
tinguished from the general model uncertainty discussed 
in section [TTl which would result when using different 
models. In order to make a quantitative estimate of this 
additional uncertainty, the spread in extracted values of 
A J (and thus F^r) was investigated by assuming specific 
forms of the interfering background missing in the VGL 
model. 

An effective upper limit for F^^ is obtained by assuming 
that the background yields a constant, negative, contri- 
bution to (Jl- For each value of Q^, this background and 
the value of Aj were fit together to the data, assuming 
that the background is constant with t. The fitted con- 
tribution of the background was found to drop strongly 
with increasing Q^. A second possibility is to assume, 
besides the VGL amplitude, a t-independent interfering 
background amplitude, fitting for every the magni- 
tude and phase of the latter, together with the value of 
Al. Although the fitting uncertainties are very large, 
the results suggest an interfering amplitude whose mag- 
nitude decreases monotonically with increasing Q^. In 
this case, the interference between the background am- 
plitude and the VGL amplitude, which depends on their 
relative phase, does not necessarily result in a net nega- 
tive cross section contribution to ctl- 

The estimated model uncertainty is determined from 
the spread of the A^ values and their uncertainties at each 
Q^, obtained with these two choices of background. To 
keep the number of degrees of freedom the same in both 
cases, the background was fixed to the value giving the 
best x^, and A J and its uncertainty were then determined 
in a one-parameter fit of the VGL model plus background 
to the data. Since the statistical uncertainties of the data 
are already taken into account in the first given uncer- 
tainty in TablelH the contribution of the statistical uncer- 
tainties of the data were quadratically removed from the 
A^ uncertainties given by the fit. The model uncertain- 
ties at each are then taken as the range plus corrected 
fitting uncertainty given by these two methods, relative 
to the value of A^ determined from the extrapolation to 
tmin- This procedure was applied to all data except those 
of Ref. [l^ , yielding the model uncertainties listed as the 
second (asymmetric) uncertainty in Table IH No model 
uncertainty was calculated for the Q^=0.35 GeV^ data 
from DESY because the t-range spanned by those data 
(only 0.03 GeV^) was too small for this procedure to be 
reliably applied. 

For the W = 1.95 GeV data, the model uncertainty in 
the extracted F^r value drops from about 20% at Q^=0.60 
GeV^ to about 5% at 1.60 GeV^. To be consistent, the 
same procedure was applied to the W = 2.22 GeV data, 
which yielded model uncertainties of about 5% at both 
Q^—1.60 and 2.45 GeV^. These rapidly dropping uncer- 
tainties with increasing refiect the smaller discrepancy 
of the VGL calculation with the t-dependence of the data 
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(8) 



where r'^ono = 0.431 fm^ is their best-fit squared pion 
charge radius. Fig. [5] indicates a systematic departure 
of the data from the monopole curve above « 1.5 
GeV^. This departure may have implications for theoret- 
ical approaches that assume the validity of the monopole 
parameterization over a wide range of . 

To illustrate the departure from the monopole curve, 
as well as to provide an empirical fit that describes the 
data over the measured range, we also show in FiglH] 
a fit which includes a small dipole component. 



F- 



fit 



85%F„ 



lb%Fdip, 



(9) 



where 



F, 



FIG. 6; (Color online) Q^F^ data from this work, compared 
to previously published data. The solid Brauel et al. [13 ] 
point has been reanalyzed as discussed in the text. The outer 
error bars for the JLab data and the reanalyzed Brauel et 
al. data include all experimental and model uncertainties, 
added in quadrature, while the inner error bars reflect the 
experimental uncertainties only. Also shown is the monopole 
fit by Amendolia et al. [^ as well as a 85% monopole-|-15% 
dipole fit to our data. 



at larger values of and W. These findings are at least 
compatible with the idea that resonance contributions, 
which presumably have a form factor that drops rapidly 
with Q"^, are responsible. They also suggest that our 
Ftt extraction methods are robust, when the background 
contribution is small, as appears to be the case at the 
higher value of W. 



dip 
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(10) 



and r^p 



0.411 fm^. This dipole parameterization 
has nearly the same for the elastic scattering data as 
the monopole curve shown [§], but it drops much more 
rapidly with . The combined monopole plus dipole fit 
is consistent with our intermediate data, while main- 
taining the quality of fit to the elastic scattering data. 
Since a monopole parameterization does not converge to 
the pQCD asymptotic limit (Eqn. [J), it is expected to fail 
at some point. Similarly, we should expect this empiri- 
cal monopole-l-dipole parameterization to show its limi- 
tations when additional high data become available 



IV. COMPARISON WITH MODEL 
CALCULATIONS 



D. Discussion and Comparison with Empirical Fits 

The form factors extracted from the Fpi-1 and Fpi- 
2 data with the use of the VGL model are shown in 
Fig. [SI along with the reanalyzed Q^=0.70 GeV^ data 
of Ref . [l3l , the elastic scattering measurements of Ref. 
[3, and the Q2=0.35 GeV^ data of Ref. [H. The Cor- 
nell data of Refs. [l^, [HI [13] are not included because, 
as discussed in section ITTB| they have large unknown sys- 
tematic uncertainties. The excellent agreement between 
the Q^=1.6 GeV^ form factor values obtained from our 
W — 1.95, 2.22 GeV data, despite their significantly dif- 
ferent trnin and W values, indicates that the model uncer- 
tainties from the use of the VGL model seem to be under 
control, at least in this Q^-range. Also shown is a more 
recently obtained value at (5^=2.15 GeV^ [i^], which was 
also extracted with the use of the VGL model. 

The solid line shown in Fig. [S] is the monopole fit 
obtained by Amendolia et al. [9| from their elastic scat- 



The pion form factor can be calculated relatively easily 
in a large number of theoretical approaches which help 
advance of our knowledge of hadronic structure. In this 
sense, F^^ plays a role similar to that of the positron- 
ium atom in QED. Here, we compare our extracted -FV 
values to a variety of calculations, selected to provide a 
representative sample of the approaches used. 

A. Perturbative QCD 

The most firmly grounded approach for the calcula- 
tion of Ftt is that of pQCD. The large behavior of 
the pion form factor has already been given in Eqn. [T] 
By making use of model-independent dimensional ar- 
guments, the infinitely-large behavior of the pion's 
quark wave function (distribution amplitude, or DA) is 
identified as 

(j)^{x,Q'^ ^ (x) ^ 6fT,x{l - x) (11) 
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whose normalization is fixed from the tt"* 



decay 

constant. Eqn. [5] follows from this expression. 

Neither of these equations is expected to describe the 
pion form factor in the kinematic regime of our data, and 
so much effort has been expended to extend the calcula- 
tion of Ft; to experimentally accessible Q^. In this case, 
the pion DA, (j)T^{x,Q'^), must be determined at finite 

. Additional effects, such as quark transverse momen- 
tum and Sudakov suppression (essentially a suppression 
of large quark-quark separation configurations in elas- 
tic scattering processes) must be taken into account. A 
number of authors 0, [4^ [3, S^l have performed 
leading-twist next-to-leading order (NLO) analyses of F^^ 
at finite Q^. The hard contributions to ^FV expand as a 
leading order part of order and an NLO part of order 

Bakulev, Passek-Kumericki, Schroers and Stefanis [i^l 
have investigated the dependence of the form of the DA 
on the form factors, using data from a variety of exper- 
iments. These were the 7r77 transition form factor data 
from CLEO [48| and CELLO as well as our F^ data. 
Their results are insensitive to the shape of the DA near 
X — 1/2, while its behavior at x = 0, 1 is decisive. The 
resulting hard contribution to the pion form factor is only 
slightly larger than that calculated with the asymptotic 
DA in all considered schemes. The result of their study, 
shown as F^"-^"^ in Fig. [71 is far below our data. The drop 
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FIG. 7: (Color online) The data of Fig. |6] are com- 
pared with a hard LO-I-NLO contribution by Bakulev, Passek- 
Kumericki, Schroers and Stefanis [ij based on an analysis of 
the pion-photon transition form factor data from CLEO "S] 
and CELLO [l^]. A soft component, estimated from a local 
quark-hadron duality model, is added to bring the calculation 
into agreement with the experimental data. The band around 
the sum reflects nonperturbative uncertainties from nonlocal 
QCD sum rules and renormalization scheme and scale ambi- 
guities at the NLO level. 



at low is due to their choice of infrared renormaliza- 
tion, which is not necessarily shared by other calcula- 
tions. To bring the calculation into agreement with the 
experimental data, a soft component must also be added. 
The treatment of the soft contribution to the pion DA is 
model-dependent. The authors estimate this soft contri- 
bution using a local quark-hadron duality model. This 
soft estimate, along with the sum of the hard and soft 
contributions, are also shown in Fig. [T] 

The interplay at intermediate between the hard and 
soft components can be non-triviaL as demonstrated by 
Braun, Khodjamirian, and Maul Q, using a light-cone 
sum rule approach. Their results support a pion DA 
that is close to the asymptotic expression, but they find 
that strong cancellations between soft terms and hard 
terms of higher twist lead to the paradoxical conclusion 
that the nonperturbative effects in the pion form factor 
can be small, and the soft contributions large, simultane- 
ously. Because of complications such as these, different 
theoretical viewpoints on whether the higher-twist mech- 
anisms dominate F,^ until very large momentum transfer, 
or not, remain. 



B. Lattice QCD 

Unlike QCD-based models, in which confinement must 
be explicitly added. Lattice QCD allows calculation from 
first principles. However, while lattice QCD is based on 
the QCD Lagrangian, it involves a number of approxima- 
tions. Errors are introduced because space and time are 
crudely discretized on the lattice. This error is controlled 
by the use of improved lattice QCD actions. To allow a 
more rapidly converging action, and hence reduce CPU 
usage, the pion mass used is significantly larger than the 
physical pion mass. Chiral extrapolation errors are in- 
troduced when the lattice results, determined with large 
pion mass, are extrapolated to physical values. Finally, 
quenching errors are introduced when disconnected quark 
loops are neglected. 

The first lattice simulations of F^^ were done in the 
1980's [sol [sil . [5^ . These pioneering works were pri- 
marily a proof of principle of the lattice technique, and 
were restricted to Q^< 1 GeV^. These results are con- 
sistent with the low experimental data, within the 
large statistical uncertainties of these pioneering calcu- 
lations. Spurred by advances in CPU power and lattice 
techniques, as well as the availability of new experimental 

data, a number of groups '53', '5^, [Hi, [Hi, [13, [Hi, [H, [13, 

have returned to the calculation of i^^ on the lattice. 
Of these, we compare our data to the recent unquenched 
simulations of Brommel, et al. [BOj. They performed 
simulations for a wide range of pion masses and lattice 
spacings, so that both the chiral and continuum limits 
could be studied. However, the lowest pion mass used 
in the simulations was 400 MeV, so the chiral extrapola- 
tion is significant. The authors fitted the Q^-dependence 
of each lattice configuration with a monopole form for 
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FIG. 8: (Color online) The data of Fig. [6] are compared 
with the lattice QCD result of Ref. [60l ] and the dispersion 
relation result of Ref. [gJ]. The lattice QCD band denotes 
the statistical and chiral extrapolation uncertainties in the fit 
monopole mass to the simulated data. The dispersion relation 
uncertainty band reflects different assumptions on the distri- 
butions of zeroes in the complex s-plane, with the 'no zeroes' 
curve lying close to the 'minimum JV' limit. 



the pion form factor and determined the corresponding 
monopole mass. They then extrapolated these masses to 
the one corresponding to the physical pion mass to obtain 
a chiral monopole mass value of 0.727 ± 0.016 GeV. The 
(monopole) form factor calculated with that mass (in- 
cluding its uncertainty) is indicated by the shaded band 
in Fig. [HI cut off at the highest point of the lat- 
tice simulation. This result begins to trend away from 
the > 1.5 GeV^ experimental data. It remains to be 
seen how these results would be affected by our Sec. IIII Dl 
comments on the applicability of the monopole parame- 
terization in this range. 



C. Dispersion Relation with QCD Constraint 

Dispersion relations are based on constraints posed by 
causality and analyticity, and relate the timelike and 
spacelike domains of the pion form factor on the com- 
plex plane. In principle the technique is exact, but our 
incomplete knowledge of the scattering amplitudes over 
the whole complex plane, and in particular the incom- 
plete understanding of the contribution of all of the poles 
in the timelike region, creates uncertainties. Authors 
address these uncertainties by imposing additional con- 
straints, such as the role of higher timelike resonances like 
the p'" , or chiral perturbation constraints near the space- 
like threshold, or that F^^ must approach its expected 
asymptotic value at very large |6^ . [63l . [s^ . [65l [6^ . 



We compare the data to the dispersion relation anal- 
ysis of B.V. Geshkenbein [U in Fig. [51 The displayed 
uncertainty band is obtained by assuming different dis- 
tributions of zeroes in the complex s-plane. This results 
in a band that grows with Q^, with the 'no zeroes' curve 
lying nearly at the lower end of the band. Our highest 
data lie above the 'no zeroes' curve, but below the 
'improved maximum ^V' limit. 



D. QCD Sum Rules 

The QCD sum rule approach is designed to interpolate 
between the perturbative and non-perturbative sectors 
using dispersion relation methods in combination with 
the operator-product expansion. While the practical im- 
plementation of this approach cannot claim to be rig- 
orously derived from QCD, its intuitive value is that it 
provides a brid ge b etween the low- and high-energy prop- 
erties of QCD [63|. A number of authors have applied 
this technique with good success to the pion form factor 
d, i, [H [el, [73. In the calculation of Radyushkin 6^, 
QCD sum rules were used to give a local quark-hadron 
duality estimate of the soft wave function 



Fr-^*(g2) = 1 



l + 6so/Q^ 

(l+4so/Q2)3/2^ 



(12) 



where the duality interval, sq, which within the QCD 
sum rule approach is determined by the magnitude of the 
quark and gluon condensates, was taken as 47r^/^2 w 0.7 
GeV^. This soft calculation, shown in Fig. [9l under- 
estimates the data by about 25%. For the hard contribu- 
tion, a simple model based on the interpolation between 
the behavior near = (related by the Ward iden- 
tity to the 0{as) term of the 2-point correlator) and the 
asymptotic behavior was used 



1 



IT (l + Q2/2so)' 



(13) 



The sum, F^°^* + F^'^^'^, is in excellent agreement with 
the data. 

More recently, Braguta, Lucha and Melikhov Q have 
replaced the simple ansatz leading to Eqn. I12l with an ex- 
pression including explicit corrections up to 0{as). Since 
the higher-order corrections needed to apply these re- 
sults with authority to the intermediate region are 
beyond the capacity of their two-loop calculation, there 
is a model dependence in their numerical result, which is 
reflected in the two different curves for sq — 0.65 GeV^ 



and So 



1+Q.(Q=)/ 



— shown in Fig. [Sj 



E. Bethe-Salpeter Equation 

The Bethe-Salpeter equation (BSE) is the conventional 
formalism for the treatment of relativistic bound states. 
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In this formahsm, a meson is described by a covariant 
wavefunction, which depends on the four momenta of its 
constituent quarks. Although formally correct, compli- 
cations arise as the interplay between different configu- 
rations, such as q-q and q-q-g, are implicitly buried in 
the potential and scattering amplitudes used in analyz- 
ing hadronic processes, and as a result, these potentials 
and scattering amplitudes are nearly intractable. The 
light-front Bethe-Salpeter model is a means to handle 
this problem by breaking the BSE into separate hard 
and soft components. A variety of models incorporating 
a confining potential which dominates at low Q^, and a 
QCD-based interaction which dominates at high Q^, are 
given in Refs. 

Another approach is to use the Dyson-Schwinger equa- 
tion (DSE) to obtain dressed quark propagators which 
may be used in the solution of the BSE. The Dyson- 
Schwinger approach to nonperturbative QCD has many 
advantages. It is consistent with quark and gluon confine- 
ment, it automatically generates dynamical chiral sym- 
metry breaking, and the solution is Poincare invariant. 
In the work of Maris, Tandy, and Roberts, the meson 
Bethe-Salpeter amplitudes and quark-photon vertex are 
obtained as solutions of the homogeneous and inhomo- 
geneous BSE, and the dressed quark propagators are ob- 
tained from the quark DSE. The model parameters are 
fixed by requiring /tt and m^r to be in good agreement 




Q2 (GeV^) 



FIG. 9: (Color online) The data of Fig. [6] are compared 
with the QCD Sum Rules calculations of Refs. [6^ and the 
Bethe-Salpeter equation model utilizing dressed quark prop- 
agators via the Dyson-Schwinger equation of Ref. [fl] [long 
dashed]. For the calculation of Ref. [6^, three curves are 
shown: [dotted] F^'"''', [short-dashed] F^"^*, and [solid] the 
sum F^-''' -I- K'"''*. For the calculation of Ref. ^, two dot- 

dashed curves are shown; [lower] so = T+cTJEPY/tv ' I^PP^^l 
so = 0.65 GeV^ 



with the data [Tg] and then and F^^ are predicted with 
no further adjustment of parameters [ttI. [78|. Their cal- 
culation is shown in Fig. [9l It is in excellent agreement 
with our data up to Q^=1.60 GeV^. To extend the va- 
lidity of the model to higher Q^, a more complete de- 
scription that takes meson loop corrections into account 
self-consistently is required [78| . 



F. Local Quark-Hadron Duality 

Quark-hadron duality relations link the hadronic struc- 
ture information contained in exclusive form factors and 
inclusive structure functions by making strong assump- 
tions of locality [79(. While local quark-hadron duality 
is an expected consequence of QCD at asymptotically 
large momenta, it is not at all clear how well it could 
work at finite fsOj. And if it does, it may be due 
to accidental cancellations of higher twist effects. Never- 
theless, it is worthwhile to compare predictions based on 
quark-hadron duality with the measured data, especially 
since duality is expected to work better at higher Q^, in 
contrast to many other approaches. 

The approximate relationship between the pion elastic 
form factor and the pion structure function F2 — 
was found by Moffat and Snell (sH . 

[FAQ^)Y ^ j'^J^^ F^{u;)dLu, (14) 

where lo = 1/x, and the upper limit of integration is 
chosen to select the elastic contribution to the inclusive 
structure function. In appl ying this formula use is made 
of the Drell-Yan-West [S^, l83| relation, which is based 
on a field-theoretic parton model that predates QCD. It 
predicts that if the asymptotic behavior of a form factor 
is (l/Q^)", the corresponding structure function should 
behave as (1 — xY'^~^ as x — *■ 1. This leads to the pre- 
diction Ff{x ^ {l-x). 

The existence of Drell-Yan data allows a quantita- 
tive test of Eon. [141 using only phenomenological input. 
Calculations [8J, [S^ based on the leading-order analysis 
of F2 data by Ref. [s^, and the next-to- leading order 
analysis of Ref. [s^j, are shown in Fig. [TOl In both 
cases, the magnitude of the prediction is dependent 
on the value chosen for the inelastic cutoff Umax (and 
corresponding Wmax) in Eqn. 1141 Local duality is ex- 
pected not to work at lower . This is reflected in the 
poor description of the < 1 GeV^ form factor data. 
However, above > 2 GeV^, the next-to- leading order 
analysis is consistent with our data. 



G. Constituent-Quark Model 

There are many calculations using a va riety of 
constituent-quark models [H, [H, S [Ml, [ol, [H, [H, [H, 
[qgI . [qtI . [9^ . The differences in approach typically involve 
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FIG. 10: (Color online) The data of Fig. [6] are com- 
pared with the local quark-hadron duality analysis of W. Mel- 
nitchouk [s^ . [ss'l , and the constituent quark model calcula- 
tions of Refs. [92, 93]. For the duality calculation, two curves 
are shown: [short-dashed] leading-order analysis of Ref. [s^ ]. 
[long-dashed] next-to-leading order analysis of Ref. [Sal . For 
the quark model calculations by Cardarelli et al. [92], two 
curves are shown: [upper dot-dashed] point-like quarks, [lower 
dot-dashed] quarks with a monopole form factor. 



differences in the treatment of the quark wave functions, 
or the inclusion of relativistic effects. Fig. [TUl shows the 
result of calculations by Cardarelli et al. 92] and by 
Hwang [93[. Both are relativistic quark models on the 
light front. Ref. [o^l uses the effective qq Hamlitonian 
of (99I , which contains a one-gluon-exchange term and a 
linear confining term, and which describes a large set of 
meson spectroscopic data. Use of this interaction results 
in large high-momentum components, and -FV is strongly 
overpredicted (upper dot-dashed curve in Fig. [10]) . This 
can be cured in a way that is consistent with the no- 
tion of a constituent quark, by assuming a form factor 
for the latter. Taking a monopole form for the latter 
and adjusting the mass parameter so that the measured 
pion charge radius is reproduced, results in the lower dot- 
dashed curve shown. 

The model of Ref. [93] allows a consistent and fully 
relativistic treatment of quark spins and center-of-mass 
motion to be carried out. A power-law wave function is 
used, whose parameters are determined from experimen- 
tal data on the charged pion decay constant, the neutral 
pion two-photon decay width, and the charged pion elec- 
tromagnetic radius. The charge and transition form fac- 
tors of the charged pion and the branching ratios of all 
observed decay modes of the neutral pion are then pre- 
dicted. The calculation is in very good agreement with 
our Ftt data. 

Li and Riska [lOOl ] asked if the empirical F.^ data ex- 



clude the presence of a significant sea-quark configuration 
in the charged pion. They performed a constituent quark 
model calculation which was extended to include explicit 
sea-quark components in the pion wave function. They 
found that these sea-quark contributions grew with in- 
creasing Q^, because they allowed the momentum trans- 
fer to be shared by a greater number of constituents, and 
so were less-suppressed at high than configurations 
which involved only a qq pair. Although their analysis 
was mo del- dependent, they found that our data allowed 
an approximate 20 ±20% sea-quark component, with the 
data point at Q^—2A5 GeV^ providing the greatest con- 
straint. 



H. Holographic QCD 

A recent th eore tical development is the AdS / CFT cor- 
respondence |lOl| between weakly-coupled string states 
defined on a 5-dimensional anti-de Sitter space-time 
(AdSs) and a strongly-coupled conformal field theory 
(CFT) in physical space-time. The goal of holographic 
QCD models is to find a weakly-coupled theory for which 
the dual strongly-coupled theory is as close to QCD as 
possible, and so allow analytic solutions of hadronic prop- 
erties in the non-perturbative regime to be performed. 
In these models, confinement is simulated by imposing 
boundary conditions on the extra fifth dimension z 
In the "hard- wall" variant, confinement is modeled by 
a hard cutoff at a finite value z = zq = 1/Aqcd- This 
has the advantage of simplicity but produces the unphys- 
ical Regge trajectory ~ n^. The "soft- wall" variant 
replaces the hard-wall boundary with an oscillator-type 
potential, and produces the more phenomenologically re- 
alistic Regge behavior ~ n. 

Several authors hav e applied holo g raph ic models to 
the pion form factor ^M, EqI ^M, [Iq3- Complica- 
tions arise when one introduces spontaneous and explicit 
chiral symmetry breaking e ffects into the soft- wall holo- 
graphic QCD model. Refs. [lOl llOSi ] take different ap- 
proaches to this problem. Grigoryan and Radyushkin 
[lOSi ] consider only the hard- wall variant, and then es- 
timate a soft- wall correction from their previous vector 
meson study |l06j . They conclude that a full analytic 
calculation would likely follow the data only in the 
< 1 GeV^ region, while overshooting it abo ve ^ 2 
GeV^. The calculations by Kwee and Lebed [104 llOTj ] 
are numerical. Both the hard-wall and the soft-wall cal- 
culations predict charge radii that are too small, espe- 
cially for the soft-wall case (see Fig. [TT|) . By allowing 
the parameters of the soft-wall model (originally fixed 
by TTip, m,r, and /,r) to vary, they find it is possible to 
describe at either high or low Q^, but not both. 
Issues in ongoing discussions 

[iolliQ3 on the AdS/CFT 
approach include the applicability of this model to the 
larger region where partonic degrees of freedom be- 
come appreciable, and the treatment of chiral symmetry 
breaking. 
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FIG. 11: (Color online) The data of Fig. |6]compared with 
the holographic QCD model calculations by H.J. Kwee and 
R.F. Lebed [l04l ]. The curves are: [solid] "hard- wall" and 
[dot-dash] "soft- wall", both with parameters fit to m,r, rup 
and and [dash] "soft-wall" with a — 262 MeV to improve 
the fit to F„ at higher but destroying the agreement with 
the other observables. 



V. SUMMARY AND OUTLOOK 

Values for the charged pion form factor, have 
been extracted for Q^— 0.60-2.45 GeV^ from the longi- 
tudinal cross sections crL{t) for the ^li{e, e'Tr'^)n reaction 
recently measured at JLab. F^^ values were also extracted 
from older experimental data acquired at DESY. The 
Cornell data are not included in this analysis because 
these ctl were not obtained in a true L/T-separation, but 
instead by subtracting a certain assumption for ax. In 
addition, the higher data have excessively large values 

of imin- 

The form factor extraction requires the use of a model 
incorporating both the 7r+ production mechanism as well 
as the effect of the nucleon. Several approaches to extract 
from the data, including the Chew-Low extrapola- 
tion method, various types of Born Term Models, and 
newer models utilizing Regge trajectories and effective 
Lagrangians, were reviewed. By using specially gener- 
ated test data, it was found that extrapolating to the 
pole at t = m^, as is done in the Chew-Low method, 
cannot be used in practice, because there is no way to 
determine the order of the polynomial to use for the ex- 
trapolation, and because even small uncertainties in the 
measured cross sections lead to a large uncertainty in F-^ . 

From the models available for determining F^^ from 
the measured values of ctl, the VGL Regge model [l^j 
was chosen, since it contains no ad hoc parameters, and 
its validity has been well established over a wide kine- 
matic range in t and W for both electroproduction and 



photoproduction data. The VGL model gives a rather 
good description of both the t and the W dependence 
of the JLab data at values of Q^> 1 GeV^, but espe- 
cially at Q2 ^ 0.60 GeV^ the fall-off of the data with 
—t is steeper than that of the model. In the cases where 
the VGL model described well the t-dependence of the ctl 
data, the value of -FV was determined by fitting the model 
to the data. Otherwise, the value of -FV was determined 
by extrapolating the fit of the model to t = tmin ■ An ad- 
ditional 'model uncertainty' has been estimated by using 
two different assumptions for an interfering background 
that could be responsible for this discrepancy between 
the data and VGL model. The fact that the discrepancy, 
and hence the model uncertainty, is very small at higher 
values of and W suggests that effects from nucleon 
resonances play a role in the data at lower and W . 

It is stressed that the cross sections are the actual 
observables measured by the experiment, and that the 
extracted values of -FV are inherently dependent on the 
model used to extract them. The development of addi- 
tional models for the ^H(e, e'Tr"*")/! reaction would allow 
further exploration of the model dependence of the ex- 
traction of from the same cross section data. On the 
experimental front, proposed measurements [4l| after the 
completion of the JLab upgrade are expected to better 
establish the validity of any used model by investigating, 
for example, the M^-dependence of the results. 

The results for FV, extracted from our data and from 
the DESY data with the use of the VGL model, are pre- 
sented together with their experimental and model un- 
certainties. Above ~ 1.5 GeV^, these data are sys- 
tematically below the monopole parameterization based 
on the empirical pion charge radius. The data are 
also compared to a selection of calculations, including 
those based on pQCD, Lattice QCD, Dispersion Rela- 
tions, QCD Sum Rules, Bethe-Salpeter Equation, Local 
Quark-Hadron Duality, Constituent-Quark Model, and 
Holographic QCD. There has been tremendous progress 
in the theory of hadronic structure physics in the past 
decade, as evident by the many new approaches under 
development. However at present, the intermediate 
regime remains a significant challenge. Several different 
approaches concur that up to at least Q^=2.5 GeV^, the 
data are far above the estimated 'hard' (perturbative) 
contribution, and that 'soft' (non-perturbative) contri- 
butions likely dominate in this region. Data expected to 
be taken 4l| after the completion of the JLab upgrade, 
up to at least = 6.0 GeV^, are expected to indicate 
whether the higher-twist mechanisms dominate _FV until 
very large momentum transfer, or not. 
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